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We investigate two types of nonlinear interaction between collinear femtosecond laser pulses with 
power slightly above the critical for self-focusing P cr . In the first case we study energy exchange 
between filaments. The model describes this process through degenerate four-photon parametric 
mixing (FPPM) scheme and requests initial phase difference between the waves. When there are 
no initial phase difference between the pulses, the FPPM process does not work. In this case it is 
obtained the second type of interaction as merging between two, three or four filaments in a single 
filament with higher power. It is found that in the second case the interflow between the filaments 
has potential of interaction due to cross-phase modulation (CPM). 
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I. INTRODUCTION 

When two, three, four or higher number of filaments propagate in parallel and close trajectories, interflow and 
observation of extreme events are reported GLSIH. These mergers, with appearing of a strong filament, were called 
Rogue events during the filamentation process. The coupling between two filaments was experimentally observed for 
first time in [l]. The conditions for optimal coupling were obtained later from Kosareva et al. Q and the optimum is 
for filaments with power near to critical for self-focusing and small diameter of the spot ro < 3 mm. As it was reported 
recently in Jj, depending on the relative phase and the incidence angle, the filament can experience fusion, repulsion, 
energy redistribution and spirial motion. The experiments in xenon gas were performed with several gigawatts input 
peak power [5j] to obtain parallel filament strings with numbers N < 12. The laser pulse breaks up from spatially 
homogeneous beam profile into several highly localized filament strings each with pulse power slightly above the critical 
for self-focusing P cr [3, § • This is the reason to look for nonlinear optical mechanisms leading to exchange of energy 
or mergers during the process of multifilament propagation. The three dimensional localization appears similar to the 
soliton interaction in one-dimensional system as optical fibers, and based on clamping effects due to CPM [MU and 
FPPM [12] . The nonlinear interaction process in fibers strongly depends on the initial phase difference between the 
pulses. 

In this paper we propose a nonlinear vector model, where in details is investigated the role of CPM and degenerate 
FPPM processes in respect to the relative movements of laser filaments. We investigate numerically the interaction 
between optical pulses in the cases when: 1) the initial phase difference between pulses is not equal to zero A tp 0 
and 2) the initial optical pulses admit equal phases A ip = 0. Thus, by properly selected initial conditions, we 
take into account the FPPM process as addition to the CPM influence. The proposed nonlinear vector model is 
investigated numerically on the base of the split-step Fourier method. We introduce by the moment formalism 
nonlinear acceleration and potentials between the weight centrums of the pulses. 


II. THEORY: NONLINEAR POLARIZATION AND BASIC SYSTEM OF EQUATIONS 

As it was pointed in TcMTpt . the filamentation process can be described more correctly by using the generalized 
nonlinear operator 


P nl = n 2 (-E • Ej E, (1) 

which includes additional processes associated with third harmonic generation. We substitute into the nonlinear 
operator m two-component electrical vector E = (E x , E y , 0) at one carrying frequency uiq 
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where A x = A x (x,y,z,t),A y = A y (x,y, z,t) are the amplitude functions and ko is the carrying wave number of the 
laser source. 

The nonlinear polarization m generates the following components 
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where fi 2 = §ri 2 . The operator 77,2 (E • Ej E generalizes the case of Marker and Terhune’s operator, and includes to 
the self-action terms, CPM terms, FPPM terms and also additional terms associated with Third-Harmonic Generation 
(THG). 

The initial laser pulses (to > 50/s) possess a relatively narrow-band spectrum (A k z <C ko) (A k z is the spectral 
pulse width). During the filamentation process the initial self-focusing broadens significantly the pulse spectrum. The 
broad-band spectrum (A k z ~ ko) is one of the basic characteristics of the stable filament. The dynamics of broad¬ 
band pulses can be presented properly within different non-paraxial models such as UPPE 0 Il4 | or non-paraxial 
envelope equations [15|. Another standard restriction in the filamentation theory is the use of one-component scalar 
approximation of the electrical field E. This approximation though, is in contradiction with recent experimental 
results, where rotation of the polarization vector is observed 0, H3|- For this reason, in the present paper we use 
the non-paraxial vector model up to second order of dispersion, in which the nonlinear effects are described by the 
nonlinear polarization components ([3]). The system of non-paraxial equations of the amplitude functions A x ,A y of 
the two-component electrical field (0 has the form 


- 2 i 


ko dA x A A P+1 (d 2 A x n d‘ 2 A x \ _ „d 2 A x 


Vgr dt 


= A I An — 


Vgr V ^ Vsr dtdz J 


p- 


dz 2 


+klfi2 


2 (^£c "b Ay) A x exp (2iko (z ( v p h v gr ) t)) + MA X | + ^ \ A V \ J A x + ^A x A y 


= a lAu -t±l(++- 2 v , t +Sl) - ++ 


V gr dt 


Vgr \ Ot 2 


d 2 Ay 

U9r dtdz J '' dz 2 


+k < Q?i 2 - (A 2 X + A 2 ) A v exp (2 ik 0 (z — ( v p h — v gr ) t)) + MAyl 2 + -|A X | 2 j A y + -A*A 2 
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where v gr and v p h are the group and phase velocities correspondingly, /? = kov 2 r k" and k" is the group velocity 
dispersion. 

This model describes the ionization-free filamentation regime, where the pulse intensities are close to the critical 
one for self-focusing. The first nonlinear term in (U) corresponds to coherent GHz generation [l5j . The system 
(|4|) is written in Galilean frame (z 1 = z — vt\t' =t). In all coordinate systems - laboratory, moving in time, and 
Galilean, the group velocity adds an additional phase (carrier-envelope phase) in the third harmonic terms and 
transforms them to GHz ones. This can be seen directly for the system ([I]) written in Galilean frame, which determines 
the choice of coordinates. The last nonlinear term in (U) describes degenerate four-photon parametric mixing. To 
satisfy the Manley-Rowe relations of the truncated equations with a generalized nonlinear polarization of the type 
P nl = n 2 ^E ■ Ej E, some restrictions on the components of the electrical field are imposed. The condition is simple 

- the possible initial components A x and A y should be complex-conjugated fields. The conservation laws give us 
additional information on the behavior of the vector amplitude function: only components of the vector amplitude 
field A = ( A x , A v , 0), which present rotation of the vector A in the plane ( x , y), satisfy the MR conditions. That is why 
in our numerical experiments, as well as in our analytical investigations, we will use complex-conjugated components 
only. 






















3 


The system of equations (©) written in dimensionless form becomes 
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where x = x/ro, y = y/r o, z = z/r$ are the dimensionless coordinates, ro is the pulse waist, zo = v gr to is the spatial 
pulse length, a = kozo, 5 = r^/zo, 'y = ^- 2 11 2 /2 is the nonlinear coefficient and A v„,i = {v p h — v gr ) /v gr is the 

normalized group-phase velocity difference. 


III. NUMERICAL SIMULATIONS 

In the experiments on multi-filamentation two basic trends are observed. The first one is that the number of 
filaments is reduced significantly as a function of the distance [IH. The second trend is observed recently in USM 
as mergers between two, tree or four filaments in one (Rogue) wave. We think that the both processes are connected 
and they are results of different types of nonlinear interaction. Therefore we investigated in details interaction of two 
filaments. By control of the initial phase difference between the pulses it is possible to include or exclude the process 
of FPPM in the nonlinear interaction. When the initial phase difference between the pulses is not equal to zero, 
the process of FPPM starts to work and an intensive exchange of energy is observed 19]. When the initial phase 
difference of the pulses is equal to zero the process of FPPM practically does not work and the nonlinear interaction 
is due the CPM process. 

The numerical results are presented for initial conditions: 120/s Gaussian bullets with waist and spatial length 
r o = zo = 72 nm and power slightly above P cr . In this case a = 907r, 5=1, A v n i = 0.00023 and 7 € 1.5 — 3. The 



FIG. 1: (la) Energy exchange between two collinear filaments Ai and A 2 with power slightly above the critical P cr (7 = 1.5). 
The pulses are separated initially at distance 2a = 3.4 and the evolution is governed by the system of equations ©. In the initial 
conditions m the phase difference ip = 7t/4 correspond to maximal energy exchange. Due to degenerated FPPM process one 
of the filaments is amplified while the other filament enters in linear mode and vanishes, (lb) Fusion between the same pulses 
when they are with equal initial phases, i. e. A ip = 0. The similar picture is seen when the FPPM terms are excluded from 
the equations © and only interaction due to CPM is investigated. With Zdiff is denoted the diffraction length Zdiff = fcor 2 . 
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FIG. 2: Interaction between three collinear filaments Ai, A 2 and A 3 governed by the system of equations lfo] l. Two of the 
pulses are with equal initial phases, while the third admits phase difference = 7t/4 in respect to others. Fusing between the 
pulses with equal phases is seen, while the third one exchanges energy by FPPM process. 


phase difference between the A x and A y components is initially ir/2 in order to satisfy the conservation laws. We 
investigate two collinear laser pulses as two vector fields Ai and A 2 at small distance a between them. Each of the 
pulses admits x and y components: Aj = Aj tX x + Aj^ v y,j = 1, 2. The initial conditions for numerical solution of the 
system of equations ([5]) have the form 
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( 6 ) 


where A x and A y are composed of the x- and y-components of the two optical pulses propagating along different 
parallel trajectories. The phase difference between the A x and A y components is initially ir/2 in order to satisfy the 
conservation laws, while the phase difference between the pulses is denoted by Atp. By varying the phase difference 
Atp we include (and exclude, when A <p = 0) the FPPM process. The interaction of optical pulses A\ and A2 for 
7 = 1.5, Av = 1.5, 2a = 3.4 and A ip = 7 r /4 is shown on Fig. la. It is observed that the amplified pulse self-focuses 
and gets enough power to continue its propagation, while the other pulse gives out energy, enters into linear mode and 
vanishes. In this way the number of filaments can be reduced by non-linear parametric processes in \^ media. In the 
following numerical experiment (Fig. lb) we exclude the FPPM process by using initial phase difference Atp = 0. To 
verify this result we also exclude the parametric step from the the system of equations ([5]) and increase the intensity by 
factor 1/4 to keep on the critical power. The pulses start to attract each other without energy exchange and as result 
a merging and self-focusing are observed. In the both numerical experiments (with A <p = 0 or when the parametric 
step is excluded from the program) the results are similar - there is no energy exchange and the fusing between the 
filaments is clearly seen. Similar potential type of interaction by CPM was reported in optical fibers [nun. in 
the following section of this paper we calculate the nonlinear acceleration and potential between the weight centers 
of optical pulses. In the general case of few optical filaments usually there are random phase differences between 
the waves. On Fig. 2 the interaction between three pulses is presented. Two of them are with equal initial phases, 
while the third one admits phase difference A<p = 7 r /4 in respect to others. Similar dependance on the initial phase 
difference is observed: fusion between the pulses with equal phases, while the third one exchanges energy by FPPM 
process. 
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IV. MOMENT FORMALISM AND POTENTIALS 

To obtain analytical expressions of the influence of CPM on the relative moving of optical pulses we exclude the 
FPPM process and GHz generation from the system of equations ©• The basic system in this case is transformed 
to (3 + 1)1? of Manakov type 
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The scalar case of nonlinear interaction is investigated in (l2j . In this paper we will preform similar analysis applied 
to collinear laser pulses presented as vector fields A\ and A 2 . We decompose as in the previous section the vectors 
Ai and A 2 in (x, y) plane 


Aj Aj^ x x + Aj y y, j — 1,2. 


Thus, the components A x and A y in (0 become 


Ax - A\ X T A 2 ,X 1 Ay - A\ y T A 2 y . 

Let us introduce the integral of energy of A x and A y 

Pj = J J J | Aj (x, y, z, t) | 2 dU = const ; j = x,y 

where dU = dxdydz and also the integrals of center of weight in x direction of A x and A y are 
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o 2 . 

Only the second derivative ^3 in the system 0 is non-commutative operator in regard to x translation. The other 
differential operators commute with x and therefore the velocity in x direction of the center of weight can be written 
as 




A* (x,y,z,t) 
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The acceleration in x direction of the center of weight can be expressed by the following convolution integral 
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(13) 


where k ^ j. In the similar way we obtain the expressions of the accelerations in y and z directions. The total 
acceleration of the components A x and A y can be presented by the following vector sums 


v gr k 0 h 2 

3pi 


^Xj(Ax,t)^ x + (Yj(Ay,t)^ y + (Zj{Az,t)jz = 
[\Aj(x + Ax,y + Ay,z + Az,t)\ 2 V\A k (x,y,z,t)\ 2 ] dU\ j = x,y, 


(14) 


where with V = d/dx + d/dy + d/dz is denoted the gradient operator of the scalar field \A k (x,y,z,t) \ 2 and the 
expression under the integral is 
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|Aj {x + A x,y + A y,z + A z,t) | 2 V|A fe ( x,y,z,t ) | 2 = \Aj (x + Ax,y,z,t) 
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Here we investigate the simplest case of two spherically-symmetric pulses, located at arbitrary distance Ax in x 
direction. Then, since Ay = Az = 0, the convolution integrals in y and z planes are equal to zero. Therefore we 
calculate the acceleration in x direction (IT51) only. Substituting the decomposition © in m we obtain 


(a(Ax,t))? = Ci 


I A X1 (x + Ax, y, z, t) | 2 — \A V2 (x, y, z, t ) | 2 + 


dx 
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(17) 


where by (a(Ax, t))£ and (a(Ax, t))^ are denoted the accelerations of the pulses (not of the components) with 
condition (d(A x,t))^ + (a(Ax, t))^ = 0 and Ci = p = PxPy . In the case of spherically-symmetric 

functions and circular polarization (A Xi = A yi ) the acceleration of the center weights can be presented in spherical 
coordinates 


{a(Ar,t)) Xi = 2Ci 
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where Ar is the distance between the centers of weight of the pulses. Since the acceleration depends on Ar (fT51) we 
can introduce the nonlinear potential 


/»A r 

V (Ar, t) = V (0, t) — (di(Ar,t)) (dAr). (19) 

Jo 

Let us suppose that the pulses do not change their shape and spectrum during the propagation process - as it can 
be seen on Fig. 16) this assumption is correct, when the pulses are at a sufficient distance from each other. At close 
distances the acceleration and potential depend significantly on time. We use trial functions with two shapes 1) 
Gaussian profile (as in the numerical experiments above) 


A = Ai = A 2 = A 0 exp - 


x 2 + y 2 + z 2 


= A 0 exp - 


and 2) Lorentz profile (such form have the filaments in the faraway zone of propagation) 


A — A\ — A 2 


2A 0 
1 + r 2 


( 20 ) 


( 21 ) 


On Fig. 3 the graphics of the nonlinear acceleration d(A r) (1151) and potential V(Ar) between the centers of 
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a) b) 




FIG. 3: Graphics of the nonlinear acceleration fi(Ar) (1181) and potential V (Ar) (1191) between the centers of weight in cases of 
Gaussian pulses (|20l) for normalized constant 2 -k 1 C\A% = 1. 

a) b) 

a(Ar) F(Ar) 




FIG. 4: Graphics of the nonlinear acceleration fi(Ar) (1181) and potential V (Ar) (1191) between the centers of weight in cases of 
Lorentz pulses m • The potential of Lorentz pulses is approximately twice wider than the potential of Gaussian ones. 


weight in the cases of Gaussian pulses GP for normalized constant 2tt 2 CiAq = 1 are presented. The same quantities 
for Lorentz pulses Oil) are plotted on Fig. 4. The formulae of the exact solutions of the convolution integral (fl8l) and 
the potentials (fTTFl) for the both cases are 


(d'i Gaus (Ar)) = —2\/27 tA r (3 + A?’ 2 ) exp (—Ar/2); y Gaus (Ar) = 2\[2 tx (5 + Ar 2 ) exp (—Ar/2). (22) 


. T . . 32-7rAr (24 + 22Ar 2 + Ar 4 ) r . 167rAr (28 + 15Ar 2 + Ar 4 ) 

d^“(Ar) =- \ a 4 - L. V L °^* (Ar) =- \/ A 2 . 3 - 1 • (23) 

(4 + Ar 2 ) (4 + Ar 2 ) 

It is important to be pointed that the potential of Lorentz pulses is approximately twice wider than the potential 
of Gaussian pulses. As a result, the Lorentz type filaments can interact at twice longer distance than the standard 
Gaussian type filaments. In the general case, the acceleration and the potential are not stationary and as it can be 
seen from the expressions m and (1191) depend in addition on the time. That is why the spatial forms and the 
spectrums of the pulses at short distances are modulated significantly. The numerical experiments demonstrate, that 
if the pulses are separated along x direction, the forms and the k x spectrums of the both pulses become asymmetric. 

















V. CONCLUSIONS 


We have developed a vector model to describe the processes of reduction of number of the filaments as well as 
the observation of mergers and Rogue events during multi-filament propagation. It is known, that in air P = P cr 
corresponds to intensity of the laser field of the order of I ~ 10 12 W/cm 2 . The main role at these intensities play 
the nonlinear effects. The results from the numerical analysis give confidence for claiming, that the investigated 
above processes are result of nonlinear interactions due to FPPM and CPM mechanisms. Finally, by using the 
method of moments, the merging between spherically-symmetric, circular polarized filaments is presented as potential 
interaction. 
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VII. LIST OF FIGURE CAPTIONS 

Fig.l (la) Energy exchange between two collinear filaments A\ and A 2 with power slightly above the critical P cr 
(7 = 1.5). The pulses are separated initially at distance 2 a = 3.4 and the evolution is governed by the system of 
equations ([5]). In the initial conditions © the phase difference p = 7 t/4 correspond to maximal energy exchange. 
Due to degenerated FPPM process one of the filaments is amplified while the other filament enters in linear mode 
and vanishes, (lb) Fusion between the same pulses when they are with equal initial phases, i. e. Aip = 0. The 
similar picture is seen when the FPPM terms are excluded from the equations d5j and only interaction due to CPM 
is investigated. With Zdiff is denoted the diffraction length Zdiff = k^r g. 

Fig.2 Interaction between three collinear filaments Ai, A 2 and A 3 governed by the system of equations (f5|). Two 
of the pulses are with equal initial phases, while the third admits phase difference A ip = 7 t /4 in respect to others. 
Fusing between the pulses with equal phases is seen, while the third one exchanges energy by FPPM process. 

Fig. 3 Graphics of the nonlinear acceleration a(Ar) (IT81) and potential V(Ar) (flT)l) between the centers of weight 
in cases of Gaussian pulses (I2U1) for normalized constant 27 t 2 C'iAq = 1 . 

Fig. 4 Graphics of the nonlinear acceleration a(Ar) lfT81) and potential V(Ar) (flT)l) between the centers of weight 
in cases of Lorentz pulses ED. The potential of Lorentz pulses is approximately twice wider than the potential of 
Gaussian ones. 


